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We study the approach, initiated by Marinari et al., to the static inter-quarlc potential based 
on Polyakov lines of finite temporal extent, evaluated in Coulomb gauge. We show that, at small 
spatial separations, the potential can be understood as being between two separately gauge invariant 
colour charges. At larger separations Gribov copies obstruct the non-perturbative identification of 
individually gauge invariant colour states. We demonstrate, for the first time, how gauge invariance 
can be maintained quite generally by averaging over Gribov copies. This allows us to extend the 
analysis of the Polyakov lines and the corresponding, gauge invariant quark-antiquark state to all 
distance scales. Using large scale lattice simulations, we show that this interpolating state possesses 
a good overlap with the ground state in the quark-antiquark sector and yields the full static inter- 
quark potential at all distances. A visual representation of the Gribov copies on the lattice is also 
presented. 
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I. INTRODUCTION 



The inter-quark potential is a principal tool in lattice 
studies of quark confinement [H. This potential can be 
understood as composed of anti-screening and screening 
effects. In perturbation theory the anti-screening struc- 
tures are paradigm effects of non-Abelian gauge theo- 
ries and their dominance under hes asymptotic freedom. 
Screening effects lower the energy but, in perturbative 
calculations, are seen to be small compared to anti- 
screening. Non-perturbatively, however, screening due to 
light quarks starts to dominate at sufficiently large dis- 
tances implyingstring breaking. Hence, as was stressed 
for example in [3| , confinement with a linearly rising po- 
tential at large inter-quark distances is a feature only of 
Yang-Mills theory with no light quarks. 

Although the potential itself must be gauge invariant, 
it is initially unclear how to construct it since the heavy 
fermions are not gauge invariant. On the lattice one mea- 
sures the vacuum expectation value of the gauge invariant 
Wilson loop and its large-time behaviour [Ij], 

WrT^e,qy{-TV{r)} , (1) 

where r and T denote the spatial and temporal extent 
of the loop. Wilson's formulation is based on connecting 
the fermions by a path ordered exponential to produce a 
single gauge invariant object. However, this is not unique 
and in this paper we will study other ways of constructing 
gauge invariant probes of the potential. 

To clarify our notation we denote states constructed 
from the heavy fermionic fields of the theory by \qx<iy)- 
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FIG. 1: The static inter-quark potential from unsmeared Wil- 
son loops. 



These are not physical states as they are not gauge in- 
variant. Gauge invariant quark-antiquark states will be 
denoted by \QxQy)- 

As discussed above, the most familiar way to ensure 
gauge invariance is to link the fermions by a Wilson line. 
Then the Wilson loop (H]) may be found from the transi- 
tion amplitude 

WrT^ {QyQx\e'"^\QM , r=\x-y\, (2) 

which describes the creation of a static quark-antiquark 
pair at t = and its subsequent annihilation at time 
t = T, see FIG. ffl 

Inserting a complete set of energy eigenstates, ((Sj) be- 
comes at large times 

WrT - \{QyQx\Oxy)\^ex_p{~TV{r)} , (3) 

where \Oxy) signifies the lowest energy state with a heavy 
quark and antiquark fixed at x and y. For practical com- 
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putational purposes it is of course important that the 
matrix element in front is sufficiently large or, put dif- 
ferently, that the state \QxQy) has a substantial overlap 
with the true quark-antiquark ground state. Gauge in- 
variance alone does not ensure this and, in particular, 
Wilson loops do not have such a strong overlap. This 
is remedied on the lattice by the technique of 'smear- 
ing' [3, 4, 5]. 

Although Wilson loops are gauge invariant and hence 
potentially physical they cannot be interpreted as am- 
plitudes connecting separately physical, single quark and 
antiquark states. This is symbolically displayed in FIG.[T] 
by the lines connecting the external sources (small cir- 
cles). However, in the perturbative sector physically dis- 
tinguishable constituent structures are expected [6]. To 
incorporate this constituent picture into the description 
of the potential we would require a gauge invariant quark 
state \Qx) built out of the heavy, but gauge dependent, 
fermionic state \qx)- 

Such an approach was, in fact, initiated by Dirac 
many years ago in QED and during the last decade or so 
this idea has been developed into a fully fledged alterna- 
tive approach to charges in gauge theories [a, H, H, ■ 
The process of constructing a gauge invariant quark from 
the matter states \qx) is called 'dressing' and details of 
the construction will be summarised in the next section. 
We will see that the dressing for a static charge takes 
its simplest form in Coulomb gauge. In that gauge, the 
fermionic state \qx) itself captures the dominant gluonic 
contribution to the perturbative quark constituent of our 
quark-antiquark configuration. As such, the state \qxq_y) 
evaluated in Coulomb gauge has a large overlap with the 



0, 



at short dis- 



true quark-antiquark ground state, 
tances. 

This observation allows us to develop an alternative to 
the unsmeared (or 'thin') Wilson loop description of the 
quark-antiquark system. We can now make the gauge 
invariant identification that for heavy static quarks. 



\QxQy) 



\<lxqy) ■ 



(4) 



This is gauge invariant as long as each orbit of the 
gauge group has one unique representative satisfying the 
Coulomb gauge condition. Thus the expectation value 
of the Hamiltonian in the state @ should give an alter- 
native approach to the inter-quark potential, see FIG. [H 

An immediate question, which we will address in this 
paper, is to what extent does the good perturbative 
overlap found in Coulomb gauge extend to the non- 
perturbative regime? If it does, then, as well as produc- 
ing an alternative to smearing, we can start to address 
issues related to the details and range of validity of the 
constituent quark picture. 

In this paper we will test this approach to the inter- 
quark potential. That is, we will use Coulomb gauge 
fermions to measure the potential. On a lattice this 
means that we first take a Coulomb dressed QQ state 
on a time slice separated by a distance r and let it evolve 
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FIG. 2; The static inter-quark potential from dressed quarks. 



during a time interval T, see FIG. [2J This yields a cor- 
relator of two separate finite length Polyakov lines with 
link variables in Coulomb gauge. The correlation of fi- 
nite length Polyakov lines was first studied on the lattice 
by Marinari et al. in They considered SU(3) gauge 
theory and obtained an upper limit for the string ten- 
sion which overestimated the true value by roughly 50%. 
It was recently confirmed (for an SU(2) gauge theory) 
that finite length Polyakov line correlators indeed yield 
the correct value for the full string tension T^. One of 
the goals of this paper is to apply this approach not just 
to the string tension but to the full inter-quark poten- 
tial. Using both perturbation theory and lattice gauge 
simulations, we will find that the potential between two 
Coulomb dressed quarks and that obtained from Wilson 
loops are in good agreement. 

An intriguing aspect of the present investigation is the 
obvious physical fact that the constituent picture has a 
limited range of validity. Thus any interpretation of the 
inter-quark potential as that arising between two physi- 
cal quarks must break down. It has been argued in refer- 
ences d, [l0| that this is consistent with the observation 
that the gauge invariance of the dressed charges breaks 
down non-perturbatively due to Gribov copies p^ . . 
In this regime the potential defined by the finite length 
Polyakov lines can only be interpreted as arising from the 
state which corresponds to one gauge invariant, but over- 
all colourless, object. However, it is not a priori obvious 
that the residual gauge invariance of the states dH) due 
to copies will let us extract a potential, let alone the con- 
fining one. Therefore, another aim of the work reported 
here is to study how Coulomb gauge Gribov copies affect 
the inter-quark potential. One might wonder whether 
summing over such copies could mean that the potential 
vanishes (as the construction is not gauge invariant with 
respect to transformations between Gribov copies). In 
fact we shall see that we can extract the potential in this 
way. 

The structure of this paper is as follows. First we shall 
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discuss what is known about the perturbative inter-quark 
potential. This will include a discussion of the relevant 
ingredients of the dressing approach to static charges and 
a summary of their use in identifying the screening and 
anti-screening structures found in the potential. This will 
allow us to clarify the significance of Coulomb gauge for 
the description of static quarks. Then, in section 3, the 
non-perturbative regime will be investigated on the lat- 
tice. Results from the use of unsmeared Wilson loops 
will be compared with those that arise from finite length 
Polyakov lines in Coulomb gauge. In section 4 we will 
then investigate in more detail the Gribov copies that 
arise in this approach and present arguments and evi- 
dence for why the inter-quark potential is insensitive to 
such copies. We will end with some comments and con- 
clusions. 



II. THE PERTURBATIVE POTENTIAL 

A. Dressing approach 

In order to be physical charges must be invariant under 
gauge transformations which map A — + ^^A and \qx) 
^^x|<7a;}- The generic form for a charged single- fermion 
state is thus given by^ 

\Q^) =h^{A\\q^) , (5) 

where the field dependent dressing, must trans- 

form as, 

h^^A\=h^\A\^i, (6) 

for IQcc) to be a gauge singlet. Note that the dressing 
h\A\ can be viewed as a field dependent gauge transfor- 
mation. The invariance of ^ from ([S]) is a minimal re- 
quirement and there are many ways to construct such 
invariant states. The simplest manner to write down 
a gauge invariant state involving a fermion- antifermion 
pair at different points x and y is to connect them by 
a straight Wilson line. This is, however, not physically 
attractive and, as the need for smearing on the lattice 
shows, is also not the best thing to do in practice. 

At a more fundamental level, the problem with such a 
description may already be seen in QED [l^ where the 
potential between two static fermions, separated by a dis- 
tance r — \x — y\ and connected (dressed) by a string, 
is easily calculated to be proportional to e^r(5^(0). This 
is a linearly rising potential (with a divergent coefficient) 
which would imply confinement in QED. However, upon 
quantisation the associated state is clearly infinitely ex- 
cited and thus energetically unphysical. On a discrete 



lattice the coefficient would be finite, but this would still 
correspond to a highly excited state which will decay. 
The reason why the energy is so high is that the electric 
field is only non-zero on the path of the Wilson line, e.g. 

E{x) = S^{x^)9ix3)n, (7) 

where n is the unit vector in the xa-direction and xj_ = 
{xi,X2)- The large energy of such thin strings is exactly 
the reason why smearing Wilson lines is numerically effi- 
cient in lattice calculations of the inter-quark potential: 
it improves the overlap with the ground state [jj 0) [E| • 

Instead of this we look for a physically motivated de- 
scription of the dressing of two static charges^. We de- 
mand that a dressed state is gauge invariant and static 
as an asymptotic state d, Q . This produces a dressing 
which factors into two parts: 

h = e'^e'^ . (8) 

For a static charge, x is the field dependent transforma- 
tion into Coulomb gauge. This part of the dressing we 
also call the 'minimal' dressing since the state e^^^'^'' Iq^) 
is gauge invariant on its own. Hence, its QQ analogue, 

|QxQj,) = e'>^('")e-^x('')|(7,g-j,) , (9) 

may be expected to have a large overlap with the true 
QQ ground state \Oxy) (involving the full dressing). The 
additional part of the dressing, k, is the time integral 
of the X rotated Aq component of the potential. So in 
Coulomb gauge it is just the (long) Polyakov line of Aq 
from the current time to infinity. 

It turns out to be a nontrivial task to find a closed 
expression for the transformation into Coulomb gauge 
in the non-Abelian case. It is, however, simple to obtain a 
perturbative expansion for x to high orders. Absorbing 
the coupling, g, into the vector potential, the minimal 
dressing at leading order (LO) is 

X^^ + 0{g^) . (10) 

In the Abelian limit this is the result proposed by 
Dirac Q to describe a static electron. At next-to-leading 
order (NLO) and beyond we have 

x-xi+X2 + X3 + 0(.g'), (11) 

where explicit expressions for the XnS can be found 
in [l3]. Given this solution for x^ it is then reason- 
ably straightforward to construct the additional part of 
the dressing to the same order in the coupling by using x 
to gauge rotate Aq. This perturbative construction can 
be extended to an arbitrary order, however, the above 
suffices to calculate the NLO potential. 



^ Note that in order to conform to lattice conventions, the gauge 
transformations have been renamed h etc. as compared 

to 0, @, [^, [loll . ^ The generaUsation to moving charges is discussed in detail in 
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A fermionic field dressed by ([8]) has two independently 
gauge invariant structures: /t and the minimally dressed 
fermion. Given these separately gauge invariant terms, 
it is instructive to investigate their relative contributions 
to the potential. The potential between two such dressed 
quarks is found by sandwiching the Yang-Mills Hamilto- 
nian between states of a dressed quark and antiquark 
spatially separated by a distance r. As in QED, the LO 
dressing generates the Coulombic electric field typical of 
short distances Q and the additional part of the dressing 
makes no contribution here. So to LO we find 



(12) 



Here, Eq collectively denotes contributions independent 
of the distance r. As usual, Cp denotes the quadratic 
Casiniir in the fundamental representation. The result 
corresponds to one-gluon exchange at tree level as 
will be briefly discussed in the following subsection. 

At NLO we focus on the minimally dressed quark 
where we only retain the Coulomb dressing. The higher 
order terms in the minimal dressing (jlip modify the es- 
sentially Abelian result p2|) and generate the paradigm 
anti-screening effects which underly asymptotic freedom. 
In momentum space [l6| , working in an arbitrary covari- 
ant gauge^ in d spatial dimensions, the chromo-electric 
commutators produce at NLO the correction 



Sg'^CpCAhkj f dfip 1 



{2tj:Y {k - p)' 



(13) 



where iDj^ {p) is the gauge invariant two-point function 

of the transverse components of the spatial vector 
potential and Ca is the quadratic Casimir in the adjoint 
representation. The potential from the minimal dressing 
at this order in perturbation theory thus becomes 
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487r2 V k^ 



(14) 



where /i is a renormalisation scale. The factor of 12 in 
front of the logarithm demonstrates that the minimal 
dressing description yields the dominant anti- screening 
part of the quark potential which is thus generated by 
the glue needed to make quarks gauge invariant. The 
expected factor of 11 in the true ground state (familiar 
from the beta function) requires an additional —1 term 
which is produced by screening due to gauge invariant 
glue — recall that the additional dressing, k, is gauge 
invariant, see [l3|- The sign of this smaller effect shows 
that it lowers the overall energy. 



^ Working in an arbitrary covariant gauge in this way allows us to 
check the gauge invariance of the construction. It also allows us 
to avoid the need for Dirac bracket commutators that arise in, 
for example, Coulomb gauge. 



In three dimensions [l8| the minimal dressing gener- 
alises this anti-screening/screening divide in the inter- 
quark potential despite the fact that this super- 
renormalisable lower dimensional theory has a vanishing 
beta-function. It is noteworthy that the relative weight- 
ing of the divide at NLO is almost identical in both 
three and four dimensions. The renormalisation of these 
physical descriptions of charges has been investigated at 
NNLO and a cancellation of non-local divergences was 
shown 19]. This was a highly non-trivial check of the 
construction. 

We note also that the arbitrary path dependence of QQ 
states dressed by a Wilson line may be factorised order 
by order in perturbation theory. This yields a product 
dressing consisting of a gauge invariant, but path de- 
pendent, gluonic term and a separately gauge invariant 
structure of two fermions dressed with the above minimal 
(or Coulomb) dressing [20l |. 



B. Wilson loop approach 

To make the transition from the dressing description to 
the lattice formulation, it is helpful to briefly recapitulate 
the perturbative route to the potential in terms of Wilson 
loops, Polyakov lines and Feynman graphs as initiated 
in pl| and reviewed in [2^ . A path integral definition 
of the static inter-quark potential V{r) may be given as 
follows. We consider the partition function 



Z[p] ^ j'^A exp{- [Sym + (P, AO)] } 



(15) 



in the presence of a static external qq source (at x and 
y, respectively). 



ply{z)^-gT-[5{z-x)-5{z~y)] 



(16) 



The measure in (|15p implicitly contains gauge fixing 
terms as does the action 5ym- The partition function 
Z is the exponential of the Schwingcr functional. 



=exp{-zl^[p]} , 



(17) 



which yields the static potential in the large-time limit, 
T — > oo, plus constant self-energy contributions E, 



W[p\^T[V{r) + Y] . 
If we define the untraced Polyakov line 



Pt{x) = Texp . 



-ig I dtAQ{t,x) 



(18) 



(19) 



we may rewrite the amplitude (|15p as an expectation 
value, or more properly a two-point function, given by 



Z[p] = {iy Pt{x)pI{v)) = J 'DAtrPT{x)P^{y), 



(20) 
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In contrast to the dressing approach to the potential, 
gauge invariance is not manifest here since the nonlo- 
cal operator tr PT{x)P^{y) transforms nontrivially as the 
Polyakov lines do not trace a closed loop. However, for 
some gauge choices (e.g. Landau or Coulomb gauge) one 
expects the spatial part of the gauge potential to decay 
sufficiently fast in temporal direction that one can close 
the integration contours to form a rectangular Wilson 
loop [23]. Choosing spatial and temporal extent r and 
T, respectively, the resulting Wilson loop (expectation 
value) is 



rT 



WrT = (trPexp<{ -ig (p dXf.A'' } ) ^ exp{-TV{r)} , 

(21) 

with the final expression being approached in the large 
time limit, T — > oo. From this familiar gauge invari- 
ant construction, perturbative calculations in covariant 
gauges of the potential have been performed to two 
loops [13, [11. 

To connect with the minimal dressing approach, it is 
useful to recall the calculation of the perturbative inter- 
quark potential in Coulomb gauge [2^,[23,[2l,[2§| . Taking 
into account that the Coulomb gluon propagator has a 
transverse and an instantaneous part, we have, at NLO, 



V 



(22) 

where the dotted lines represent 'unphysical' instanta- 
neous Coulomb gluons, Ag, with propagator 



^00 (fe) 



(23) 



while the 'curly' lines correspond to 'physical' transverse 
gluons Af with propagator 



k) 



(24) 



The full vertical lines represent the heavy-quark propa- 
gators or, equivalently, the untraced Polyakov lines. 

Looking at we note that only gluonic two-point 
insertions contribute. The second diagram corresponds 
to the minimal contribution to the potential (fT5|) while 
the third diagram is the smaller screening contribution 
and may be traced back to the separately gauge invariant 
dressing, k, of 

Higher n-point functions will only contribute at two- 
loop order and above in the coupling. We mention in 
passing that this order has never been worked out in 
Coulomb gauge as one is faced with a number of diffi- 
culties. First, there are severe infrared singularities [30| . 
Second, integrations over energies (p") tend to diverge 
badly. Third, there are nontrivial contributions from 



the Faddeev-Popov measure and the Christ-Lee terms 
[31I [3^ induced by the curvature of the configuration 
space. Items two and three seem to be related (33l[33.[35|. 

Working out the diagrams of ((22l) the potential in mo- 
mentum space becomes 



K.o(fe) 



(25) 



where the quantity in brackets determines the running 
coupling. Asymptotic freedom arises because the anti- 
screening or minimal-dressing contribution yields a factor 
-f 12, cf. (dH), as compared to the —1 from the screening 
vacuum polarisation term. 

One can understand the different signs at NLO by not- 
ing that 'physical' particles (like the transverse gluons) 
in a loop always entail screening by unitarity '36, 37] . At 
small distances this effect is due to virtual particles in the 
loop. At large distances, however, screening effects may 
be caused by real particles. This is the case for string 
breaking in QCD where light dynamical quark pairs pop 
out of the vacuum and combine with the heavy probe 
quarks to form heavy-light mesons. This in turn leads 
to a saturation of the potential at long distances. The 
same happens in pure gluodynamics if the heavy sources 
are in higher group representations such that they can 
be screened by dynamical gluons. 

Beyond NLO the separation into screening and anti- 
screening structures is not so straightforward and the 
only systematic approach is through the dressing decom- 
position ([H]). So, for example, at NNLO physical particles 
in loops can now contribute to anti-screening [l9|. The 
gauge transformation properties of the dressing ensures 
that, at all orders in perturbation theory, the separation 
of forces into screening and anti-screening ones is gauge 
invariant and hence physical. How this is achieved in the 
non-perturbative sector is an open question. 



III. THE NON-PERTURBATIVE REGIME 

The use of Coulomb gauge in the non-perturbative 
regime has received much attention recently: there have 
been investigations using lattice techniques |13 . ISa . [39l . 
liol . l4ll . [H, [4^, variational arguments [43, |45|. l46l. [47f. 
Schwin ger- Dyson equations [48j] and explicit analytic so- 
lutions [10]. In particular the so-called Coulomb poten- 
tial has been studied in great detail. This potential arises 
from the short Polyakov line correlator (|20p in the limit 
T ^ [il [m, HI m and must not be confused with 
the full static potential, which we address in the present 
study. In this section we shall recall the key ingredients 
of the lattice approach and then see how to extract the 
inter-quark potential from finite length Polyakov lines. 
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A. Numerical setup 

The dynamical degrees of freedom are the unitary ma- 
trices Ufj^{x) S SU(2) which are associated with the hnks 
of a iV"* cubic lattice with lattice spacing a. The partition 
function is given by 

J T)U p[U] , with p[U] = exp {/3 S[U]} , (26) 

where the Wilson coupling /3 is the only free parameter of 
the simulation. It multiplies the gauge invariant Wilson 
action, 

(27) 

which is used throughout this paper. The integration 
measure is given in terms of single-site Haar measures 
dC/p(x): 

'DU = \{dU^,{x). (28) 

The latter are invariant under left and right multiplica- 
tion by SU(2) group elements and normalised to unity, 
J dU^{x) — 1. This ensures the gauge invariance of DU 
and, therefore, of the partition function. Under a gauge 
transformation ^{x), the links change according to 

Uf^ix) '^U^{x) = n{x) U^ix) n\x + fi) . (29) 

In order to impose Coulomb gauge in a lattice simulation, 
one firstly determines a link dependent gauge transforma- 
tion, h{x), from the 'action principle', 

SR^ = J2^^^U,{t,x) max. (30) 

which, on each gauge orbit, looks for the representative 
^Ui that is closest to the unit matrix. Note that the sum 
over times, t, guarantees that this is done term by term 
on each time-slice. 

Lattice configurations in Coulomb gauge are then 
made up from the gauge transformed links ^U^{t,x). 
The search for a maximum of the gauge fixing func- 
tional S'fix is performed with a standard iteration and 
over-relaxation procedure. For reasonable lattice sizes, 
present day algorithms will typically find a local (rather 
than the global) maximum of the gauge fixing action ([50]) 
on each orbit. 

Starting with the same 'seed' configuration C/^(a;), we 
might end up in several, almost degenerate, maxima of 
([50)1 implying that the lattice definition of the Coulomb 
gauge configuration is ambiguous. This is the lattice ver- 
sion of the continuum Gribov ambiguity discussed in the 
introduction. Two configurations ^^U^{x) and ^^Ufi{x) 
with both hi{x) and h2{x) satisfying ([50)1 are hence called 
Gribov copies of each other. Since we are primarily in- 
terested in the effects induced by Gribov copies, we do 



not employ sophisticated tools such as simulated anneal- 
ing or Fourier acceleration which are invoked (with little 
success it must be said) to find the global maximum of 

dsni). 

B. Confining potential from unsmeared Wilson 
loops 

The standard way to calculate the potential of a static 
quark-antiquark pair is to investigate the transition am- 
plitude (21), which we now write as 

Wir,T) = {QyQx\e-"^\Q^Qy) , (31) 

where, as before, H denotes the Yang-Mills Hamiltonian, 
T a Euclidean time interval and r = \x—y \ the separation 
between quark and antiquark. Recall that the latter are 
connected by a parallel transporter along a straight line. 
We will evaluate W{r, T) via lattice Monte Carlo simu- 
lations employing rectangular Wilson loops (see FIG. [J). 
As discussed in the introduction, there are sophisticated 
lattice techniques 0, 0, Q available for enhancing the 
overlap, {QyQx\Qxy) , with the quark-antiquark ground 
state by using (spatially) smeared Wilson loops. How- 
ever, since we are not only interested in the static poten- 
tial per se, but also in the overlap of different interpo- 
lating states with the true ground state, we refrain from 
adopting these techniques and content ourselves with a 
thorough comparison of unsmeared Wilson loops, where 
there is only an overall gauge invariant quark-antiquark 
state, with individually dressed Coulombic states as pic- 
torially summarised in FIG.s[T]and[21 

Let us first discuss the evaluation of PT|) . Rewriting 
the amplitude as 

M^(r,r) =exp{-r?;(r,T)}, (32) 

the inter-quark potential V{r) emerges as the limit 

V{r) = lim v{r,T) . (33) 

T — >co 

The large T limit ensures that excited states which 
contribute to the matrix element l|3ip are suppressed. In 
this limit v{r,T) becomes independent of T. FIG. [3] (up- 
per panel) shows — lnW{r,T) = Tv{r,T) as a function 
of T for several values of r. Deviations from the linear 
behaviour are clearly visible at small values for T. The 
curves for larger values of r seem to be more affected by 
excited states. 

In practice, we fit the logarithm of ([5^ . 

^liiW{r,T) = b{r) + V{r)T (34) 

to a linear function in T (we call this a T-fit). Only data 
points with T > Tth, for some threshold time Tth, are 
taken into account. Using small values for the thresh- 
old time, say Tth = 2a, the x^/dof becomes unaccept- 
ably large. This indicates large deviations from linear 
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FIG. 3: The dependence of unsmeared Wilson loops on T; 
16* lattice, /? = 2.4. 



behaviour and implies that the signal is contaminated 
by contributions from excited states. We have therefore 
made a more conservative choice for the threshold, 



Tth = 4 a 



(16'' lattice, /? ^ 2.4) (35) 



and included only T-fits satisfying 
Xr/dof < 2 . 
The lower panel of FIG. [3] shows 
Vsir,T) = [- In W{r,T) 



(36) 



6(r)]/T = u(r,r) - &(r)/T 
(37) 

as a function of T. It is reassuring to note that, for 
T > 4a, the contributions from excited states are within 
the statistical error. The T-fit described above also pro- 
vides us with the desired static inter-quark potential 
l^(r). Generically, the quark-antiquark pair is aligned 
along the main crystallographic direction, e.g. {100}, of 
the lattice. In order to check for artifacts arising due 
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FIG. 4: The static potential extracted from unsmeared Wil- 
son loops as a function of the quark-antiquark distance r using 
Tth = 4a (upper panel) and Tth = 3a (lower panel) as thresh- 
old for the T-fit; 16* lattice, /3 = 2.4. 



to rotational symmetry breaking, the spatial sides of the 
Wilson loop were placed along the diagonal axes {110} 
and {111} as well. Our final result for the unsmeared 
loop and Tth = 4a is shown in the left panel of FIG. ID 
To increase the number of statistically significant data 
points for V(r) we also relaxed the condition ([55)1 by us- 
ing the threshold Tth = 3a, see FIG. SI lower panel. Also 
included in this graph is the solid line from the fit to the 
Tth = 4a data. We find that Tth = 3a data approach 
this line reasonably well. Note that the contribution of 
the excited states in the Tth — 3a data yields a slight 
overestimate of the string tension. 

Finally, the potential y(r) is fitted to the universal 
function 



y{r) = h err , 

r 



(38) 



with offset Vq , Coulombic coefficient a and string tension 
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TABLE I: Fit of the static potential to the function (f38l) for 
a 16'' lattice and for f3 = 2.4. 
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a 


cra^ 


Xv/doi 


Tth = 3a 


0.487(1) 


0.204(1) 


0.0908(4) 


13.0 


Tth = 4a 


0.501(3) 


0.212(2) 


0.0847(8) 


4.7 



a to be determined (we call this a T^-fit). Our findings 
are summarised in Table ID Not unexpectedly, we find 
that caution is required when the static potential is ex- 
tracted from unsmeared Wilson loops: There is a signif- 
icant contribution from excited states the suppression of 
which requires a sufficiently large Euclidean time extent 
T. For the present settings employing a 16'* lattice and 
f3 — 2.4, we find a 7% drop in the string tension from 
aa"^ w 0.0908 to aa'^ « 0.0847 if we impose the more 
rigorous bound ([55)1 for the T-fit. Note also that calcu- 
lations which use smeared Wilson loops for an overlap 
enhancement find smaller values for the string tension, 
see e.g. [i^. This indicates that, even for the case of the 
more conservative bound ([55]) . there are probably still 
some contributions from excited states. 



C. Confining potential from Coulomb dressed 
quarks 

In what follows we will employ lattice gauge config- 
urations which have been transformed to the Coulomb 
gauge along the lines described in subsection IIII Al Our 
working hypothesis here is that we can extract the phys- 
ical intcr-quark potential despite the presence of Gribov 
copies. The authors of [111, found that the large dis- 
tance part of the potential, i.e., the string tension, is 
not affected by Gribov copies. However, no physical ex- 
planation of this observation was given. Although the 
copies obstruct the interpretation of the resulting poten- 
tial as that arising between separately invariant single 
quark sources ffl we shall demonstrate that they do not 
affect the overlap with the ground state. 

The overlap of a Coulomb dressed, single heavy quark 
state at i = with the analogous state at t = T is given 



where the links can, in fact, be taken in any gauge. Thus, 
C(r, T), and with it the potential U{r) from ((4T|) . are not 
affected by the trivial residual gauge freedom ([1^ . 

For a quark trial state which is not completely invari- 
ant because of Gribov copies, one might naively expect 
that the average over gauge copies causes the correlator 



by a finite length Polyakov line (evaluated in Coulomb 
gauge), 

T 

^PT{x)^l['^i^^Uoix,t), (39) 

t=0 

where ^U^{x) denotes the link after gauge fixing. Accord- 
ing to ([20]) the potential of two static Coulomb dressed 
quarks can then be extracted from the correlator of two 
Polyakov lines of finite temporal extent, 

C(r,r) = (tr''PT(a3)''4(y))- (40) 
This correlator was studied in [l^ with an emphasis on 
the Coulomb (minimal anti-screening) potential obtained 
from the limit T ^ 0. In this paper, we will concen- 
trate on the full static potential (large T behaviour) and 
the impact of Gribov copies. For sufRciently large val- 
ues of T the excited states are suppressed and we expect 
C(r, T) to be dominated by the true QQ ground state 
Oxy)- Hence, defining quantities analogous to p2|) and 

(El, 

C(r,r) = exp{-ru(r,T)} , U{r) = lim w(r,r), 

(41) 

we expect U{r) to coincide with the static potential V{r) 
extracted from unsmeared Wilson loops. The purpose of 
the present subsection is to scrutinise this expectation by 
a detailed numerical investigation. 

As we will have to deal with gauge fixing ambigui- 
ties let us first address the issue of residual gauge in- 
variance. As is appropriate for a Hamiltonian formula- 
tion, the Coulomb gauge prescription ([30]) is implemented 
on time slices, t = const. Hence, we may still perform 
spatially homogeneous, purely time dependent residual 
gauge transformations, ^{t), on the gauge fixed links, 

''u^it.x) -> ''^'^''u^{t,x) = n{ty'u^{t,x) n\t), (42) 

which, in continuum language, obviously leave the (van- 
ishing) divergence of A{t^x) invariant. 

Let us now convince ourselves that the correlator 
C{r,T) of (gni) is invariant under (li^ . 



(43) 

I 

C{r,T) to vanish. Our numerical simulations, however, 
show that this does not happen: we obtain a clear signal 
for C(r, T), which is orders of magnitude larger than the 
underlying noise. We therefore apply the linear T-fit to 



I 

tr "WpT(a;) '"^''^pUv) = f^(0) Pt{x) n\T) n{T) P^iy) n\0) = tTPrix) P^iy) , 
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TABLE II: Fit of the static potential U{r) to the function 
(pHj) for a 16* lattice and for /3 = 2.4. 



FIG. 5: u{r, T) extracted from the correlation of finite length 
Polyakov lines calculated from Coulomb gauge fixed lattice 
configurations (upper panel). The static potential calculated 
from finite length Polyakov lines (lower panel). Both on a 16* 
lattice with (3 = 2.4. 



the function — lnC(r, T) subject to the rigorous bounds 
Tth = 4a , xr/dof < 2 . (44) 

The fit results allow us to extract u{r,T) to high ac- 
curacy as shown in FIG.O upper panel. The outcome is 
clearly encouraging: not only is the Gribov noise absent, 
it also appears that the overlap with the ground state is 
largely enhanced. In FIG. \5\ even the r = 7a data are 
statistically significant. This is also true for the plot of 
the static potential U{r) defined in (|4ip . Finally, we have 
again fitted the parametrisation ([38| of the potential to 
the present data (FIG. [5l lower panel). The results of 
this F-fit are listed in Table HIl For comparison, the pa- 
rameters for the fit with threshold Tth — 3a have also 
been included. The results seem to be more stable than 
those from the fits to unsmeared Wilson loop data (cf. 
Table |T| . This once again signals the larger ground state 
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Xl/doi 


Tth = 3a 


0.505(1) 


0.214(1) 


0.0827(2) 


11.0 


Tth = 4a 


0.510(2) 


0.217(1) 


0.0807(4) 


6.5 



overlap of the Coulomb dressed QQ states as compared 
to the qq states linked by an unsmeared Wilson line. 

Concluding this subsection we reemphasise the strik- 
ing fact that Gribov copies seem to have no noticeable 
effect on the inter-quark potential. As a result, possible 
contamination by excited states notwithstanding, it is 
fair to say that the potentials calculated from unsmeared 
Wilson loops and from Coulomb dressed quark states do 
agree. 



IV. ON THE IMPACT OF GRIBOV COPIES 

In the last section we have shown that the static po- 
tentials obtained from finite length, Coulomb dressed 
Polyakov lines and from unsmeared Wilson loops agree. 
Our calculations have required a substantial processing 
of raw numerical data such as the subtraction of diver- 
gent self- energies. The question therefore arises whether 
the agreement between standard and Coulomb gauge cal- 
culations also holds for the bare finite time amplitudes, 
W{r,T) and C(r,T) from ^ and (gT]), respectively. 

If this were true, we would be in a much better position 
to understand the impact of Gribov copies as this issue 
could be properly addressed before renormalisation. 



A. Gribov noise 

Let us start with a description of the actual numeri- 
cal procedure which is used to calculate the expectation 
value of an observable 0[U] in Coulomb gauge. Using 
standard Monte-Carlo techniques, a link configuration 
{[/} is generated according to the gauge invariant prob- 
ability distribution p[U], cf. Without gauge fixing 
an expectation value is calculated as 



(O) = / VUpiUpiU] 



(45) 



If we evaluate this after an arbitrary gauge transforma- 
tion, U ^ ^U, using invariance of the probability mea- 
sure T)Up[U], 

(O) = J DUp[U]OfU] = j TiUp[U] J DflOfU] , 

(46) 

we obviously obtain zero whenever the projection 
/ "DilOl^U], i.e. the group average of O, vanishes. For 
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gauge variant field combinations this is what typically 
happens. 

However, the situation will be different after gauge fix- 
ing where, as a result of the gauge fixing algorithm, any 
seed configuration will change into a Coulomb gauge con- 
figuration 



In this case, (|52|) results in the expression. 



{U} 



h[U] 



(47) 



with (yet unknown) probability P^U] characterising the 
algorithm. We emphasise that the proper way to inter- 
pret the transformation (|47|) is as a change of variables 
[3l| implying curvilinear coordinates and a metric formu- 
lation which, however, will not be needed here. 

It will be sufficient to define the expectation value of 
an operator 0\U] in Coulomb gauge according to 



m) 



(48) 



where the double bracket emphasises that 0\U] need not 
be gauge invariant. Note, however, that the dressed op- 
erator O [''[/] is, by construction, invariant up to Gribov 
copies as h always picks an orbit representative. If there 
were no Gribov noise the chosen representative would be 
unique. This uniqueness is corrupted by the unavoidable 
randomness inherent in the gauge fixing algorithm im- 
plying a distribution of representative copies. It is this 
distribution which we want to study in the following. 

To this end let us rewrite the expectation value ((48l) 
in terms of the Yang-Mills functional integral. The prob- 
ability of picking a particular copy ''[/is given by the 
gauge group average 



tt{^U\U) = / Dnp{^U\^U) , 



which is gauge invariant by fiat^ 



(49) 



(50) 



The gauge dependent quantity p{^UY''U^ in p9|) denotes 
the (gauge dependent) probability distribution for ob- 
taining the output copy when the numerical proce- 
dure was initiated with input and hence characterises 
the gauge fixing algorithm. With these prerequisites, we 
find for the distribution P in 



P[^U] = j DU t:{^U\U) p[U] . (51) 
Inserting this into (pS)) the expectation value becomes 
((O)) = / 'DUpp] f D''{/^(''{/|[/) 0[''U] . (52) 



Before we proceed further let us check that we recover 
the familiar formula under the assumption that there is 
no Gribov ambiguity. In this case, the gauge fixing algo- 
rithm would map U to the unique Coulomb gauge repre- 
sentative U implying a sharp distribution 

nCulU) ^ 6{''U\''"U) , where ho ^ HqIU] . (53) 



m) 



= J DUp[U] 0[''°U] 



(54) 



which is gauge independent by the (assumed) uniqueness 
of hop]. 

In the generic case, when Gribov copies are present, 
([55| is no longer true. Nevertheless, gauge invariance 
is maintained due to the invariance property (j50p of 
7r(''[/|C/). This allows for the definition of a gauge in- 
variant average over Gribov copies. 



0[U] = / 1)''U ttC'UIu) 0[''C/] , 



(55) 



which indeed satisfies 0[^?7] — 0[U]. Finally, upon 
plugging this into (|52p . the expectation value may be 
compactly written as 



((O)) = (O) , 



(56) 



where the brackets on the right-hand side indicate the 
standard average with the Yang-Mills probability distri- 
bution, p[U]. 

Let us apply this now to the finite length Polyakov 
loop correlator, i.e. 



0[U] = C[U] =triPTix) P^iy)) , 



(57) 



and analyse its gauge invariant average C[U] associated 
with Coulomb gauge. Since C[U] is gauge dependent 
it could nevertheless happen that this average ('Gribov 
mean'). 



C[U] = / D'^U t:{''V\U) C [''[/] 



(58) 



though gauge invariant, actually vanishes when the 
copies are averaged over. We will show shortly that this is 
not the case. Therefore, C[U] can indeed be interpreted 
as a nontrivial gauge invariant observable. 

Let us consider a specific example. For a 16"* lattice 
and f3 = 2.4 we have generated a single sample configura- 
tion to which we apply a random gauge transformation. 
The result, call it U is submitted to the gauge fixing algo- 
rithm which produces a Gribov copy configuration, ''^[7. 
With this configuration, we calculate the correlator 



^5]tr '^^Pt{x) ^^Pliy) 



(59) 



where y = x + re^. We choose the distances r = 3a 
and r = 4a for the present example. Repeating this 
procedure 200 times, leaves us with observables Ci, i — 
1 . . . 200. Their numerical values are distributed as shown 
in FIG.[6l Note that for a gauge invariant (i.e. noise free) 
field combination, such as the plaquette, one would find 
a delta function distribution according to (|53p . However, 
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n □ 



OIL 



0.046 0.048 



0.05 

c 



0.052 



0.054 



FIG. 6: Gribov noise distribution of the correlator C from 
(|57|) for one sample configuration U and distances r = 3a and 
r = 4a. 



for C\U] in (j58p we find a broad distribution with mean 
/i and standard deviation s given by: 



0.0498(1) , and s = 0.0011(1) 



(60) 



A bootstrap analysis was used for the error estimate in 
these values. As they depend significantly on the gauge 
fixing algorithm as well as the observable under consider- 
ation their interpretation in terms of physical quantities 
is somewhat limited. They do, however, provide infor- 
mation on the impact of Gribov copies since (i) /i signif- 
icantly differs from zero and (ii) the relative width s//z 
of the distribution is only of the order of 2%. This rules 
out the possibility that the average over Gribov copies 
induces a null result. 



B. Gauge invariant signals from Gribov- noisy data 

As a first step, we will investigate the relation between 
the data from finite length Polyakov lines and gauge in- 
variant Wilson loops. Again using the 16^ lattice and 
(3 = 2.4, we have generated 10 seed configurations U . For 
each seed configuration, we produced 10 Gribov copies of 
it. We have then calculated r — ia and T = Aa Wilson 
loops which are averaged over the spatial coordinates. 
The result is called Vl^[?7]. For each Gribov copy, the 
spatial average of finite length Polyakov lines was also 
obtained (denoted as C^U]). As in the previous sub- 
section, there will be a distribution of values of C 
because of the Gribov noise. 

This data is shown in FIG. El For each seed configura- 
tion (and therefore for each value M^[J7]) there is a band 
of values for C^U] produced by the Gribov copies. If we 
assume a perfect correlation of the expectation values, 
i.e. 



u ' 




0.024 0.026 0.028 0.03 0.032 0.034 0.036 0.038 

W[U] 




FIG. 7: Scatter plot for ten seed configurations U (upper 
panel). Scatter plot of u{r,T), see (|4T|, versus v{r,T), see 
(|32p . for several values of r and T (lower panel). 



the data would be symmetrically scattered around the 
solid line also shown in FIG. U\ If, on the other hand, 
the Gribov noise were uncorrelated to the underlying seed 
configuration, the scatter plot would be homogeneous. 
Our findings, shown in FIG. [7] left panel, suggest that a 
significant correlation survives the Gribov average. 

To further test this correlation between finite length 
Polyakov lines and gauge invariant Wilson loops, we com- 
pare v{r,T) from ([5^ as obtained from unsmeared Wil- 
son loops, with u{r,T) from (|1T|) . FIG. [71 lower panel, 
shows the corresponding scatter plot. We clearly observe 
a linear correlation between the two quantities: 



u(r,T) = v{r,T) + e{T) . 



(62) 



((C)) (X (M/) 



(61) 



The important observation is that e{T) does not con- 
tribute to the r dependent part of the static potential. 
Hence, we find that only the unobservable offset of the 
potential is changed if finite length Polyakov lines are 
employed rather than unsmeared Wilson lines. This ob- 
servation has far reaching consequences: Gribov copies 
will affect the quark self-energies, but their impact on 
the r dependence of the static potential is limited. 
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C. The Gribov gallery 

In order to trace out the origin of the ideal correlations 
discovered in the previous subsection, let us visualise the 
spatial dependence of the Gribov copies. This depen- 
dence can change the r dependence of the finite length 
Polyakov line correlator as follows. We note that a sim- 
ilar investigation for the case of Landau gauge was per- 



formed in [51| . Generally two different runs of the gauge 



C/^. This in turn yields 
'^'^Uf^, which are related 



fixing algorithm result in two different maxima along the 
orbit of a seed configuration, 
two Gribov copies, say '*C/^ ar 
by a residual copy gauge transformation, c. Knowing 
that the Coulomb gauge ([50)1 has a trivial residual gauge 
freedom consisting of purely time dependent gauge trans- 
formations (j42p . we need to distinguish between genuine 
copies ch and those of the form n{t)h{x) which would 
also satisfy the Coulomb gauge condition. We therefore 
define 



c{x) = n{t)c{x) 



(63) 



and use the purely time dependent degree of freedom fl(t) 
to bring c as close to unity as possible. 



tr c(x, t) — > max 



(64) 



This should suppress the contamination by trivial copies. 
In other words, a space dependent c{x) with tr c{x) /2 not 
being 1 everywhere signals the presence of a non-trivial 
Gribov copy. With these prerequisites we can now inves- 
tigate how such a copy affects the finite length Polyakov 
loop correlator. 

Defining the spatial pure gauge link. 




FIG. 8: Region of space where trc(a;)/2 resulting from (|64p 
significantly deviates from 1, both on a 16* lattice (upper 
panel) and on a 24'* lattice (lower panel) 



M{yx,t) = c\y,t) c{x,t) 



(65) 



the Polyakov line correlator associated with the copy 
^'^Un may be written as 



xy L 



- (tr M{yx,0) ''Pt{x) M\yx,T) ''P^(y)) , 

(66) 



This means that, whenever the matrices c{x) differ signif- 
icantly from unity, an additional non-local correlation is 
induced to the correlator of finite length Polyakov lines. 

FIG. [5] illustrates the space dependence of c{x) for two 
generic lattice configurations for /3 = 2.4 at a fixed time 
slice. The density of points in the graph is a direct mea- 
sure for the deviation of trc/2 from 1. Although the 
underlying seed configuration is UV noisy, trc changes 
rather smoothly throughout space. This explains why 
M [yx) is not too different from unity for moderate dis- 
tances \y — x\. Note that M{yx) strongly depends on the 
efficiency of the gauge fixing algorithm in singling out a 
unique maximum on the gauge orbit. Without adapting 



the algorithm, there will be also a substantial dependence 
of M{yx) on the volume since the algorithms generically 
become inefficient for large volumes. This fact can be 
seen from FIG. [5] where two different volumes have been 
used in conjunction with the same gauge fixing algorithm. 
For the larger volume (lower panel), the modulation in- 
duced by M[yx) is much more pronounced. 

An attractive feature of the Coulomb gauge is the fact 
that the Gribov link, M{yx,t), is purely spatial (i.e. lo- 
cated in a single time slice t) and thus does not inter- 
fere with time evolution and the transfer matrix formal- 
ism. Hence, as long as the field combination C[''[/] does 
not vanish upon averaging over Gribov copies, this av- 
erage only affects the definition of the interpolating trial 
state, \QxQy)- In other words, different gauge fixing 
algorithms will correspond to different ansdtze for the 
quark-antiquark trial states which might differ in their 
overlap with the true ground state, \Oxy)- This is some- 
what reminiscent of the situation arising in the process 
of smearing Wilson lines. 
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V. CONCLUSIONS AND DISCUSSION 

In this paper we have developed a formalism based 
on Coulomb gauge fixing which allows us to incorporate 
the ubiquitous Gribov copies in a systematic and gauge 
invariant way. This has been achieved by defining, for 
the first time, a Gribov average of gauge variant field 
combinations. 

We have applied this to a thorough lattice investi- 
gation of Polyakov lines of finite temporal extent in- 
troduced by Marinari et al. in [ll[ and further studied 
in [13, SH, 0, E^. These lines, even if traced, are not 
often studied on the lattice since their expectation value, 
without gauge fixing, is zero due to their lack of gauge 
invariance. Using our formalism, we have shown that, 
and explained why, they have a physical interpretation 
when working in Coulomb gauge. This has enabled us to 
extract the confining inter-quark potential and to anal- 
yse the Gribov copies associated with that gauge and 
their impact on the overlap with the ground state in the 
quark-antiquark sector. 

We have seen that the value of the string tension 
obtained from the correlator of finite length, dressed 
Polyakov lines lies below that obtained from our simu- 
lations with unsmeared Wilson loops (compare Tables 
1 and 2). This indicates that the Coulomb dressed 
Polyakov lines are, as we expected, less contaminated by 
overlaps with higher energy states. Indeed, even using 
shorter lines in this approach gives a smaller string ten- 
sion than for unsmeared Wilson loops of larger temporal 
extent. 

Many topics brought to light here deserve further ex- 
ploration. In particular, Gribov copies have been demon- 
strated not to hinder the measurement of the full quark- 
antiquark potential in this construction. They merely 
contribute to the definition of the interpolating quark- 



antiquark state, and we have demonstrated in this paper 
that its overlap is even better than that of the state where 
quark and antiquark are joined by a thin flux line. It 
would be interesting to study in a more systematic way 
the contributions of excited states to the finite length 
Polyakov lines. There is also an obvious extension of 
this, namely to use other gauges and so test whether, as 
we predict, other dressings have a poorer overlap with 
the true ground state of Yang Mills theory. 

The Gribov induced non-locality seen in section 4.3 
hinders a proper definition of the constituent quark pic- 
ture. We would argue that our results support the, ex- 
pected, breakdown of a constituent picture of quarks and 
antiquarks in the confining region. The structure of the 
copies displayed in FIG. [51 their volume dependence and 
sensitivity to algorithms therefore merit further investi- 
gation. It would also be very interesting to probe the 
Gribov non-locality in the high temperature deconfine- 
ment regime where the constituent quark picture is ex- 
pected by many to make sense. In this context, it would 
also be important to study the relation between the Gri- 
bov non-locality and the confining vortices advocated in, 
e.g. [S^. If the confining vortices are removed, see [SOj . 
do the Gribov copies seen in FIG. [5] vanish with them 
thereby giving rise to a constituent quark picture? 
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